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Abstract 

Explicit relations between the isospin mixing probability, the spreading width 
t\ as of the Isobaric Analog State (IAS) and the statistical decay width T c 
of the compound nucleus at finite excitation energy, are derived by using the 
Feshbach projection formalism. The temperature dependence of the isospin 
mixing probability is discussed quantitatively for the first time by using the 
values of r|^ s and of T c calculated by means of microscopic models. It is 
shown that the mixing probability remains essentially constant up to a tem- 
perature of the order of 1 MeV and then decreases to about 1/4 of its zero 
temperature value, at higher temperature than ~ 3 MeV, due to the short 
decay time of the compound system. 
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Isospin symmetry is the first dynamical symmetry in nuclear physics proposed by W. 
Heisenberg in 1932 [I] and has played a very important role in the study of the structure of 
atomic nuclei for many years. One of the highlight of the isospin symmetry is the discovery of 
the isobaric analog state (IAS) of the parent ground state in the daughter nucleus by means 
of charge exchange reactions [[J. Recently, problems related to the isospin symmetry and 
its breaking have been discussed with new emphasis because new experiments in nuclei near 
the drip lines have become possible. These radioactive beams facilities can provide access to 
regions of large isospin deformations || . Several theoretical predictions have been made for 
the isospin mixing and large mixing probabilities in heavier nuclei near the proton drip line 
are found, up to (5-6)% in 100 Sn Possible effects of these large isospin impurity are at 
present under consideration, both by experimentalists and theoreticians, in connection with 
the width of the IAS and the Fermi f}~ decay of N=Z nuclei. 

A restoration of the isospin symmetry in highly excited compound nuclei was pointed 
out by H. Morinaga |J and D.H. Wilkinson || already 40 years ago. This problem has been 
studied experimentally by measuring the cross sections for the de-excitation of the giant 
dipole resonances (GDR) in thermal, self-conjugate nuclei. The experimental data show the 
sharp decrease of the mixing amplitude at high excitation energies above £'*=20 MeV in 28 Si 
and 60 Zn f/J as well as in 40 Ca f8j. The restoration of the isospin symmetry was theoretically 
discussed by using a simple model in Ref. ||. 

In this letter, we would like to address the problem of the restoration of the isospin 
symmetry by using a microscopic model, which was recently proposed and used at zero 



temperature in Ref. |10| . This model is based on the Feshbach projection method and 
gives an explicit relation between the spreading width of the IAS and the isospin mixing 
probability. In Ref. this relation was employed successfully to calculate, from the values 
of the isospin mixing, the spreading widths of the IAS in Sn- and I-isotopes and in 208 Bi. 
In the present work, we will extend the microscopic model to the case of excited compound 
nuclei, and study quantitatively the temperature dependence of the isospin mixing through 
its relation with the spreading width of the IAS. 



The isospin mixing amplitude is written in perturbation theory as | ] I 



« To+1 = rp 1 ciTo+1 (M;T = T + 1,T Z = T \H 1 \T = T ,T Z = T >, (1) 
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where H\ includes the isospin non-conserving interactions, i.e., the isovector part of the 
Coulomb interaction and the charge symmetry breaking as well as charge independence 
breaking interactions. The state with isospin To is the mother state (either the ground state 
or the IAS), while the states with isospin To + 1 give rise to the isospin impurity in the mother 
state by coupling to it. In the cases of medium-heavy and heavy nuclei, the isovector giant 
monopole resonances (IVGMR) give the major contribution to the isospin mixing so that 
we consider only the coupling to these states in the following. 

At finite temperature the mother state lies at some high excitation energy and its lifetime 
as well as that of the IVGMR built on top of it, are affected by statistical emission of 
particles (mainly neutrons) and 7-rays. The corresponding width is the compound width T c 
and the denominator of Eq. ([I]) should contain the sum of the imaginary terms iTj° /2 and 
iT^ 0+1 /2 |TJ|]. We simply denote in the following this sum as iT c /2 and T c will be called 



"compound width". In general, we put a tilde over the widths at finite temperature. The 
compound width will be added to the sum Ym of the intrinsic escape and spreading widths 
of the T + l monopole state. We discuss later how all these widths are estimated by using 
microscopic models. 

Let us now study the relation between the isospin mixing amplitude and the spreading 
width of the IAS. We will use the same Feshbach projection method of Ref. |10[ by extending 



it to the case of excited compound nuclei. The whole model space of nuclear configurations is 
built up with states which are eigenstates of the isospin-conserving nuclear Hamiltonian Hq 
and is divided in two subspaces P and Q. The P space contains only two states, the parent 
ground state and the IAS which have isospin T . The other nuclear states in the daughter 
nucleus which have isospin T — 1 form the space Q. The spreading width of the IAS (f \ AS ) 
is obtained by coupling with doorway states of the space Q through the isospin-breaking 
Hamiltonian as |10| 



fLs(£*) = -2/m(IAS | H 1 Q-^——QH 1 | IAS) (2) 

E - H QQ 

where E* is the energy of the IAS in the compound nucleus and Hqq is the matrix ele- 
ment of the isospin-conserving Hamiltonian H between states in Q-space, among which the 
dominant role is played by the three isospin components of the IVGMR in the daughter 
nucleus. At finite temperature, both the states in P- and Q-spaces have finite decay times, 
as discussed above. In terms of T c and Tm, Eq. (0) is rewritten as 

| (IAS \H l \q) | 2 



T l IAS = -2ImY: 



E*-E q + tT c (E*)/2 + iT M (E*)/2 



= V(T (E*) + f (E*)) 1 (IAS | Hi | q) | 2 
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where the E q are the excitation energies of the three different isospin components of the 

IVGMR in daughter nuclei, which have T'=T — 1, T and T + 1 and T C (E*) means, as 

recalled above, the sum of the two compound widths at the excitation energy corresponding 

to the nulear temperature, and at that energy plus the IAS energy E*. The isospin structure 

of the IVGMR can be treated properly by the Feshbach projection method as it has been 

done in Ref. flQfl . The right hand side of Eq. (||) results as a product of three factors, by 



assuming that for T c and Tm average values over q can be taken, and precisely 

f\ AS = (T C (E*) + f M (iT))(a To+1 ) 2 (To + 1)F(T ), (4) 

where the isospin factor F(T) is given by 

1 2T-1 1 {T-lf 1 4T 2 1 2 

F{T) - t {{ 2TTT> (AElf 1 ) 2 + T + l {AE T M f + (2T + 1)(T+1) (AEj/ 1 ) 2 ^ ' 

(5) 

In this last Equation, AEf.j = E\[ — E* is the energy difference between the three isospin 
components of the IVGMR in the daughter nucleus and the energy of the IAS whereas 
AEmw — E M +1 — E w is the energy difference between the T + l component of the IVGMR in 
the mother nucleus and the mother state \tt). This last difference enters because the isospin 
mixing amplitude a To+1 is obtained by Eq. ([!]). In Ref. ]10|] the isospin splitting was taken 



care of by using the Lane potential and we write a similar expression for nuclei at finite 
temperature, 

AEXj = E* (IVGMR) + -j-t ■ T c , (6) 

where E*(IVGMR) is the average energy of the IVGMR in the compound nucleus and 
t(T c ) is the isospin of the IVGMR (the daughter nucleus). The value E* (IVGMR) is taken 
from the systematics in nuclei at zero temperature as 

E* (IVGMR) = 170 A~ 1/3 MeV. (7) 

The value of the parameter of the Lane potential commonly used is V\ = 120 MeV. In 
general, the Brink-Axel hypothesis suggests that giant resonances built on excited states lie 
at the same energy as the corresponding resonances built on the ground state. In particular, 
in the case of the isovector giant dipole resonance it has been experimentally observed that 



the excitation energy does not vary appreciably with increasing temperature |l3fl . Thus, it 
is expected that the energy of the IVGMR also does not vary with increasing temperature. 

The isospin geometrical factor F(T Q ) is certainly important to predict a proper isotope 
dependence of the spreading width f \ AS of several nuclei in medium-mass region. For our 
present purpose, however, we may be allowed to discard the isospin dependence of the 
IVGMR energies. Then, Eq. (|j) becomes 

f \ AS - (f c (E*) + f M (E*))(a T ° +1 ) 2 (T + 1) (8) 

so that the isospin mixing probability is given by 



[a 



T +l\2 



(9) 



T + ir c (E*) + r M (E*) 

In Eq. (||), T c is expected to rise significantly with temperature. On the other hand, the val- 
ues of the spreading width f j AS and of the total width f m are not expected to be very much 
temperature dependent, because of the following reasons. The temperature dependence of 
the intrinsic escape widths of giant resonances is discussed in Ref. [13] and is found to 



be small. For the intrinsic spreading widths IN-, the authors of Ref. Jl5] predict a very mild 



temperature dependence. Consequently, from Eq. (|9|) we should expect at high temperature 
a strong decrease of the isospin mixing, or equivalently a restoration of the isospin symmetry. 
Eq. (|9|) certainly justifies the conjecture which was raised by Wilkinson || to describe the 
isospin mixing problem: if the compound nucleus decays on a time scale which is shorter 
than the time needed for a well-defined isospin state to mix with states with different isospin, 
then the isospin symmetry is partially or totally restored. Consequently, the restoration of 
the isospin symmetry can be understood on the ground of simple kinetic arguments (as 



pointed out also in Ref. ||). The same kinetic arguments are used in Ref. [16], where an 
equation with the same physical content as our Eq. (|9|) is obtained within the context of 
the study of the disappearance of the collective giant resonance strength at high excitation 
energy. 



The isospin mixing in compound nuclei was also discussed extensively in Ref. ||17|| . The 
authors derived a formula for the isospin mixing probability o? in self-conjugate nuclei as 



(« To+1 ) 2 = % (10) 

1 2 

where Y\ and Y\ are the mixing and total decay width of the states in the channel 2 which 
is equivalent to the Q-space of our model. Although they do not specify the spin and parity 
of the states in channel 2, the formula ([T0|) looks similar to our Eq. (||) except for the isospin 
factor which is necessary for non self-conjugate nuclei. It should be noticed that the two 
models give almost the same final result, although the formalisms by us and by the authors 



of [0 look at first sight very different. 

To discuss quantitatively the dependence of the isospin mixing on the excitation energy 
of the compound nucleus, in the example of 208 Pb, we have to model the temperature 
dependence of all quantities which appear in Eq. (|9|): the intrinsic widths of the IVGMR, 
T\j and fjy-, and the intrinsic spreading width of the IAS f \ AS , as well as the compound 
width f c . First, the intrinsic spreading widths f and f \ AS are parametrized by means of 
a linear temperature dependence according to the findings of fl5| , 
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f l = T l (T = 0) • (1 + cT), (11) 

where T is here the nuclear temperature in units of MeV and c is a constant which should 
be determined using the results of a microscopic calculation. In general, giant resonances 
are described microscopically as a coherent superposition of one particle-one hole (lp-lh) 
excitations. The main cause of their spreading width stems from the coupling to collective 
low-lying surface vibrations. In Ref. 1)151 , the spreading width of the GDR in 120 Sn has 



been studied and the slope parameter c has been estimated to be about 0.025 or less, by 
fitting the results of microscopic calculations which take into account the coupling of the 
giant resonance with states composed of a lp-lh pair plus a low-lying vibration |fL8| . These 
doorway states to which the giant resonance is coupled directly, are not considered to be 
sharp but distributed over an energy interval of the order of A with a Lorentzian shape. In 
this spirit, the width A is supposed to take care of the coupling of the doorway states with 
the compound nucleus states. This coupling has been shown, on the basis of rather general 
hypothesis, to be essentially constant as a function of temperature ||19|| . On the other hand, 
we may let A grow when increasing T to simulate the fact that the doorway states are more 
fragmented since the low- lying states become less collective [18]. It is found [2(J that the 



change of the spreading width of the IAS is at most 20% when A is varied from 50 keV 
to 500 keV, and the spreading width saturates for larger values of A. If this 20% change 
of the spreading width f j AS is supposed to occur within the temperature range between 
and 4-5 MeV, the slope parameter c is c = 0.04-0.05, which is not much different from the 
value mentioned above, c = 0.025. Consequently, we take the linear dependence QIID with 
c = 0.025 for fj^f and f \ AS . Then, we show below that our conclusions are not markedly 
changed even if we increase by a factor 2 our value of c. 

At zero temperature, the spreading width of the IVGMR r] w -(T = 0) in 208 Bi, at the 
excitation energy of the IAS, is found to be of the order of 250 keV Q, while that of the 



1 This calculation was already reported in Ref. [10| where the total width of the IVGMR in 208 Bi 
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IAS T\ AS (T = 0) is evaluated to be about 100 keV |20] including the contribution of the 



charge-symmetry breaking and charge- independence breaking interactions | 21| . This value 
T l IAS (T = 0) = 100 keV is in good agreement with the experimental value Hl7| , p2| . Finally, 
T M (T) is taken as a constant according to the findings of Ref. |14| and fixed at 200 keV 0. 

The compound widths T c are dominated by the most probable decay process of the 
compound nucleus which is usually neutron evaporation. We may discuss the neutron evap- 
oration in simple terms, as it is done in Ref. [^3j where a Weisskopf formula is used for the 
temperature dependence of r„, 

f n = ^TV^, (12) 
Tin 

with a radius R = 7.1 fm and a neutron separation energy B n = 7.4 MeV in the case of 
208 Pb. It is known that the neutron widths obtained through this simple formula may be 
somewhat different from more realistic neutron widths, essentially due to the schematic 
assumption of a geometric absorbtion cross section for neutrons on the nucleus. In order 
to take into account the microscopic transmission coefficient, we have also estimated the 



compound widths by means of the statistical computer code CASCADE [24]. This code 
takes into account also proton and a emission and needs as an input the excitation energy 
E* of the compound system which can be related to the nuclear temperature by the well- 
known relation E* = aT 2 for which we have taken a =10 MeV -1 . The results for the isospin 
mixing (|9|) are not so markedly sensitive to whether we use compound widths produced by 
CASCADE or by the simple formula (|i~2D. Therefore, for simplicity we will use the values 
of T c = T n obtained by using Eq. (|12|). 

We present our results for the isospin mixing probability of 208 Pb as a function of temper- 
ature in Fig. 1. All curves start, at zero temperature, from the value a 2 = 0.1, corresponding 



at the IAS energy was estimated to be about 450 keV: about a half of the total value, that is, 250 
keV, corresponds to the spreading width whereas 200 keV corresponds to the escape width. 
2 See the previous note. 
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to a value of the spreading width of the IAS, T\ AS = 100 keV. The solid line corresponds 
to the calculation described above with the slope parameter c = 0.025. The dashed line 
corresponds to choosing the parameter c equal to 0.05 both for the IVGMR and for the 
IAS. These two results show that our general conclusions are not too much affected by the 
choice of c and that the isospin restoration starts to take place at a temperature of about 
1.5 MeV. At 3 MeV temperature, the value of the isospin mixing probability has decreased 
to one fourth of the value at zero temperature due to the rapid increase of the compound 
width appearing in the denominator of Eq. (|9]). 

Experimentally, as mentioned above, the temperature dependence of the isospin mixing 
in the decay of the GDR was studied in the light nuclei 28 Si and 60 Zn [[?]. A similar substantial 
decrease of the isospin mixing was found above T=l MeV although the magnitude of the 
mixing at low temperature is much larger than in the present case because of the coupling 
to the shell model T> states, which is more important than the coupling to the IVGMR in 
light nuclei . 



In summary, we have studied the isospin mixing in compound nuclei at finite temperature 
with a microscopic model based on the Feshbach projection method. We are able to derive 
a relation which connects three physical quantities, that is, the isospin mixing amplitude, 
the spreading width of the IAS and the statistical decay width of the compound system. In 
this way, the isospin mixing probability in 208 Pb as a function of the temperature is studied 
quantitatively for the first time by using the results of microscopic calculations. It is shown 
that until relatively low excitation energies, corresponding to a temperature of the order of 
1 MeV, the mixing probability remains constant. At higher excitation energies, this mixing 
probability decreases and around 3 MeV it is reduced by about a factor 4 because of the rapid 
increase of the compound width. Thus, we have verified, by using microscopic ingredients, 
the validity of the conjecture formulated by Wilkinson 40 years ago, of the restoration of the 
isospin symmetry in highly excited systems. In connection with the present work, it would 
be of extreme importance to obtain experimental information about the width of the IAS 
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at finite temperatures. 
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Fig.l 

The temperature dependence of the isospin mixing probability a 2 in 208 Bi. This probability 
is calclulated by using Eq. (|9]), and assuming a linear temperature dependence for f and 
f\ AS according to Eq. ([11]). The slope parameters c adopted in ( |TT| ) are 0.025, and 0.05 for 
the solid and dashed curves, respectively. The value of r c is evaluated through Eq. (|T2"D. 
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